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Abstract 

We perform a phenomenological study of the scalar sector of two models that 
generate neutrino mass at the three-loop level and contain viable dark matter can¬ 
didates. Both models contain a charged singlet scalar and a larger scalar multiple! 
(triplet or quintuplet). We investigate the effect of the extra scalars on the Higgs 
mass and analyze the modifications to the triple Higgs coupling. The new scalars 
can give observable changes to the Higgs decay channel /i —)• 77 and, furthermore, 
we find that the electroweak phase transition becomes strongly first-order in large 
regions of parameter space. 
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1 Introduction 


The Standard Model (SM) of particle physics is a spectacularly successful theory that 
stands as one of the truly great scientihc achievements. Despite this success, however, the 
theory possesses a number of short-comings, suggesting it will likely require extensions 
and/or modifications in the future. The most obvious motivation for extending the SM 
is the need to incorporate gravity. However, a lack of present-day experimental guidance 
makes the pursuit of the theory of quantum gravity an incredibly challenging task. 

Additional evidence that the SM is incomplete comes from the experimental observa¬ 
tion of neutrino mixing and the need to explain the missing gravitating (i.e. dark) matter, 
required on galactic scales. These puzzles motivate the addition of new particle species to 
the SM and there is much hope that such new species will manifest in future experiments 
(in particular at Run II of the LHC). The neutrino mass and dark matter (DM) problems 
have stimulated much research and there are multiple candidate solutions one can pursue. 
The problems may have independent solutions, though it seems reasonable to ask whether 
the two puzzles could share a unified or common solution. Could the neutrino mass and 
DM problems be related? 

In 2002, Krauss, Nasri and Trodden (KNT) proposed a simple extension of the SM 
model that admits a relationship between the existence of DM and the origin of neutrino 
mass [1]. In this approach, one adds new fields to the SM, such that neutrino mass is 
generated radiatively at the three-loop level, with one of the particles propagating in the 
mass diagram being a DM candidate. The model employs two charged singlet scalars, 
Si 2 j and three generations of gauge-singlet fermions N. A Z 2 symmetry with action 
{S*^, N} —)■ {— —N} is also imposed. This ensures stability of the lightest fermion 
N, thereby giving a DM candidate, and also prevents a coupling between SM neutrinos 
and N, which would otherwise generate tree-level neutrino masses. The result is a type 
of unihed description for the origin of neutrino mass and DM, with the removal of the 
DM candidate simultaneously turning off neutrino mass0 

In recent years, a number of basic generalizations of the KNT model have appeared. 
In one such model (hereafter ‘the triplet model’), the singlet fields S 2 and N are replaced 
by SU{2)l triplets [6]. This model retains the Z 2 symmetry to ensure DM stability and 
prevent tree-level neutrino masses, and gives a viable alternative unified framework for 
the DM and neutrino mass problems. A further generalization exchanges the singlet 
fields S 2 and N for SU{2)l quintuplet helds [7]. This model (hereafter ‘the quintuplet 

"‘For recent studies of the KNT model see Refs. laiaiiia. 
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model’) does not require the symmetry in order to prevent tree-level neutrino masses, 
and is a viable theory for radiative neutrino mass, independent of DM considerations. 
Interestingly, the most-general Lagrangian for the model contains a single Z 2 -breaking 
coupling A. When taken small, A <C 1, the model gives a long-lived DM candidate, while 
turning A off completely, A ^ 0, activates a Z^ symmetry and gives absolutely stable DM. 
Thus, the quintuplet model is a viable model of radiative neutrino mass, with or without 
DM. 

Due to the presence of larger multiplets with non-trivial SU ( 2 ) charges in both the 
triplet and quintuplet models, the phenomenology of the models is rather rich. In the 
present work, we extend the analysis of Refs. ilE] and undertake a more extensive study of 
the phenomenology of both models. We investigate the effect of the triplet and quintuplet 
scalars on both the Higgs mass and the triple Higgs coupling, showing that the latter can 
experience sizable modifications. We also study the effect of the new multiplets on the 
Higgs decay channels h —)■ 77 , yZ. Our work shows that, e.g. observable changes are 
expected to B{h —)■ 77 ), with some regions of parameter space already excluded for the 
triplet model. The effect of the enlarged scalar sector on the electroweak phase transition 
is also analyzed, revealing a tendency for a strongly first-order phase transition in large 
regions of parameter space. 

Before proceeding we note that further generalizations of the KNT model are possible. 
Ref. [ 8 ] presented colored generalizations and other related three-loop models that employ 
slightly modified loop topologies. A septuplet generalization of the KNT model was 
proposed in Ref. [S]. This had the interesting feature of automatically containing an 
absolutely stable DM candidate, without requiring a new symmetry. A minimal scale- 
invariant implementation was also recently studied [lO]. More generally, a number of 
authors have studied connections between radiative neutrino mass and DM in recent 
years, see e.g. Refs. m-m- 

This paper is structured as follows. In Section |2l we outline the triplet and quintuplet 
models, describing some key features and elucidating some stability constraints on the 
scalar potentials. We study the influence of the new multiplets on the Higgs mass and 
the triple Higgs coupling in Section [31 The electroweak phase transition is considered in 
Section H] and we turn to the Higgs decay channels h —)■ 77 and h 7 Z in Section |5l 
Conclusions are presented in Section | 6 l 
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2 Three-Loop Radiative Neutrino Masses 


The SM employs the gauge symmetry Qsm = SU{3)c x SU{2)l x U{1)y- In this work, 
we consider extensions of the SM that include the charged singlet scalar S r'U ( 1 , 1 , 2 ), the 
scalar multiplet T ~ (1, 2n + 1, 2) and three generations of chiral beyond-SM fermions, 
(1, 2 ? 7 , +1, 0), where i = 1, 2, 3, labels generations and numbers in parenthesis denote 
charges under Qsm- We use the integer n = 0 , 1 , 2 , to label the distinct models. The 
case with n = 0 is the KNT model, for which all beyond-SM fields are SU{2)l singlets: 
S = Si , T = S 2 and F = N. For n = 1 (n = 2 ) the multiplets T and F are SU{2)i 
triplets (quintuplets) and one has the triplet (quintuplet) model. In all cases, the new 
multiplets are subject to a discrete symmetry with action {T, J^} ^ {~T, —-F}. This 
ensures a stable DM candidate, which should be taken as the lightest fermion = J-biviE 
Detailed analysis of the DM annihilation channels appears in Refs. PE]. 

With the aforementioned particle content, the Lagrangian contains the following terms: 

+ - V{H,S,T). (1) 

Here, La ~ (1,2,—1) are SM lepton doublets, CaR ~ (1,1, —2) are the SM charged 
lepton singlets and fa^ = —fga denote Yukawa couplings. Lepton flavors are labeled by 
lower-case Greek letters, a,(3 E {e, /i, r}. The singlet-leptons CaR couple to the exotics 
T and T through the Yukawa matrix gia, and the superscript “c” is used to denote 
charge conjugation. The Lagrangian shows that both T and F” are sequestered from SM 
neutrinos. We denote the SM Higgs doublet as H (1.2,1). 

The combination of the Yukawa terms in Eq. ([1]) and the term ~ in the 

scalar potential V{H, S, T) (discussed below) explicitly breaks lepton-number symmetry. 
The models therefore generate radiative neutrino masses, which appear at the three-loop 
level as shown in Figure [H Due to the Z 2 symmetry, the neutral components of the 
exotic fermions F" do not mix with SM neutrinos at any order in perturbation theory, 
and similarly there is no mixing between charged leptons and T. In both the triplet 
and quintuplet models, the charged scalar S can be within reach of TeV scale collider 
experiments 

^For n = 0 the scalar T is charged, while for n = 1 and n = 2 the neutral component of the scalar T 
has non-trivial couplings to the Z boson and id therefore excluded as a DM candidate by direct-detection 
constraints. 
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Figure 1: Three-loop diagram for neutrino mass. Here, S ~ (1,1, 2) and T ~ (1, 2?7,-|-l, 2) 
are beyond-SM scalars while ~ (1,2?7, -|- 1,0) is a beyond-SM fermion. The case with 
n = 0 corresponds to the KNT model [1], while n = 1 gives the triplet model [ 6 ] and 
n = 2 gives the quintuplet model [7]. 

2.1 Triplet Model 

The case with n = 1 gives the triplet model, for which Ti ~ (1, 3, 0) and T ~ (1, 3, 2) are 
SU{2)l triplets. We write the triplet helds in symmetric-matrix notation as T = Tab and 
= J'ab, where a, 6 G {1, 2} are SU{2)l indices. The multiplets contain the following 
components [B] 

7'„ = 7'+, J-u = 7-21 = T^^ = F]: = (7*Y, 

r„ = r++. r,2 = r2i = ir+. = (2) 

while the triplet mass term gives 

- (3) 

The neutral-fermion mass terms are brought to the correct sign by dehning the Majorana 
fermions Radiative corrections from SM gauge bosons lift the degen¬ 

eracy between the components of leaving as the lightest component [16], though 
for most purposes this small splitting can be neglected [ 6 ]. We work in the diagonal basis 
with = diag(Mi, M 2 , M 3 ), where Mi = M^m is the lightest triplet-fermion mass. 
According to the analysis in Ref. [ 6 ], the DM mass should lie in the range 2.35 — 2.75 
TeV, and all triplet members masses should be larger than the DM mass. Accordingly, 
the triplet scalars are beyond the reach of the LHC, though, as we shall see, they can still 
be probed indirectly. 
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The scalar potential for the triplet model has the form 


V{H, S, T) = + \\H\^ + + ^\S\^ + + |-[(T*)“'T,,]2 

+ 'f{T*r%c{T*r'^T,a + XsMSmi^ + {AstI^P + \ht\H\‘^} [(T*)“'T„,] 


( 4 ) 


Vacuum stability requires that the quantities 


A, \s, 7]1 + r/ 2 , hi + |h 2 , 


A (Aht — ^ht) 


As X^sT 

[Xht — Xht) Vi + V 2 


XsT Vl + lV2 


(5) 


are taken strictly positive, with A^ = min(A^, 0). Taking the charged scalar squared mass 
/i| and the scalar triplet squared masses as positive ensures the absence of spontaneous 
charge symmetry breaking and guarantees that (T^) = 0. The latter is is necessary to 
preserve the Z 2 symmetry and retain a stable DM candidate. 


2.2 Quintuplet Model 

The quintuplet model corresponds to n = 2, in which case one has J^i ~ (1,5,0) and 
T ~ (1,5,2) as SU{2)l quintuplets. In symmetric-matrix notation, the quintuplets are 
written as Tatcd and J^abcd, where [7] 

77+ TTO f 77+'iC 

L+ T- _ - TT _ _ _ 


77+ 770 ( 774 

77 _ 77++ 77 _ '' L 77 _ '' L 77 _ \'' R J 77 _ 77++1C 

1111 — L 1 •' 1112 — '' ~ '' ~ ’ '' 


R ) ’ 


rp + + rp-\- rjnO 

T' _ T'+++ T' _ T' _ T' _ T' _ T' — 

.11111 — J - , .11112 — - .11122 — -‘■1222 — -‘2222 — J - ■ 

V4 V 6 \/4 


( 6 ) 


Observe that and T are distinct fields with T ^ [T^)* . The explicit expansion of 
the fermion mass term gives 

1 


(7) 

where is a Majorana fermion and the other four components of combine to give two 
charged (Dirac) fermions: 


.F++= J-+++ J-++, ^+ = J-+-J-+, ^o = j-o + (j-) 


-Onc 


( 8 ) 


Without loss of generality, we again employ the basis with Mp- = diag(Mi, M 2 , M 3 ), 
where Mi = Mom- According to the analysis of Ref. [7], the DM mass is expected to lie 
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in the range 5.65 — 6.95 TeV, and all quintuplet members should have masses that exceed 
M,m§ Consequently the quintuplet scalars cannot be produced directly at the LHC. 

The full scalar potential for the quintuplet model that respects the global symmetry 
Z 2 is given by 

T) = l^r + + >^sh\S\^\H\^ + 

T |[(TT-T..]^ + I[(TT-T..,(TT^-T...] + I[(TT-T..(TT^-T.,.] 

+ {Xst\S\^ + Xhti\H\^} + XHT2{TT^'^T,UH*rHa 

+ jiS-fTabcdTefghe^‘^e^^e'^3e‘^’^ + H.c. (9) 

Vacuum stability requires that the following quantities 


A 


A, A5, ?7i + |r72 + Vi + \V2 + |h 3 , Vi +'>12 + V3 > 0 

/, , 3 , \0 Xgrp 




^HTl + IA 


^HT2) Vl 


‘HTl + |A_h-T2) 
^2 + |h3 


A‘ 

^ST Vl + \V2 + \V3 (hi + 
^n , 1 ^ , 1 _ \0 


XgT 


1 X n 


be strictly positive. Similar to the triplet case, spontaneous charge symmetry breaking 
is avoided and the neutral quintuplet remains VEV-less by taking the squared masses of 
the charged scalar and the quintuplet to be positive. This preserves the Z 2 symmetry. 

In the numerical scans performed below, we impose the above mentioned conditions 
such as vacuum stability, charge non-breaking and (T°) = 0, and also require the Higgs 
mass to be within the range reported by ATLAS and CMS, rrih = 125.09 =F 0.21 GeV 
El- We restrict our attention to the perturbativity domain, demanding that the physical 
vertices in Eqs. ([4]) and ([9]) be less than 3. For both the triplet [6] and quintuplet [7] 
models, we consider the charged singlet scalar mass between 100 GeV and 1 TeV. The 
masses for the scalar multiplet members should be larger than the dark matter mass, i.e., 
Mt > 2.35 TeV and Mt > 5.65 TeV, for the triplet and quintuplet models, respectively. 
For the numerical analysis we consider 20,000 sets of benchmark points for both the triplet 
and quintuplet models. The benchmarks reproduce the observed DM relic density while 
also achieving viable neutrino masses and avoiding lepton flavor violating constraints (see 
Refs. O E] for discussion on constraints). 

®The degeneracy between neutral and charged components of J- is again lifted by radiative corrections. 
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3 Higgs Mass and Triple Higgs Coupling 


In order to estimate the Higgs mass and the triple Higgs conpling at one-loop, it is 
necessary to properly dehne the effective potential, with the Higgs mass being its second 
derivative and the triple Higgs conpling given by the third derivative: 






d‘i 




h='L 


T 


( 11 ) 


Here h is the real part of the nentral component in the donblet, v is its VEV, and 

is the zero temperatnre one-loop Higgs effective potential. In this work we employ the 


DR scheme, for which the effective potential is given by [18 


yr=0/' 


('>) = - +!'>“+E 


A 


Ui 




647r^ 


log 


A 2 


( 12 ) 


where n* is the held mnltiplicity and A is the renormalization scale, which we take as 
the measnred value of the Higgs mass, A = 125.09 GeV HZ]. The quantities mf{h) 
are the held-dependent squared masses (presented in the appendix for both triplet and 
quintuplet models). In this class of models, h is the only scalar with a non-zero VEV, so 
all held-dependent masses can be written as (h) = _|_ a‘fh‘^/2, for constant ai. 

At tree-level, the parameter /i^ in the potential is given by /i^ = Xv^. After including 
one-loop corrections, the parameter /i^ is corrected as 


-I- 


327r2 


E 


njCKim, 


O 

mf 






(13) 


in order to ensure the one-loop VEV value remains as n = 246 GeV. The term 
represents the radiative corrections to the /i^-term, due to all helds, and is expected to be 
dominated by contributions from the new heavy helds. The Higgs mass at one-loop can 
be similarly dehned by 


ml = 2Xv^ + 


V 


327r2 


^ ruaj log 


A 2 ’ 


(14) 


where the radiative corrections (i.e., second term in flT4jl I are also expected to be domi¬ 
nated by contributions from heavy new helds. Although ml is determined by experimental 
observations, the doublet quartic coupling A can still be (very) small, relative to the SM 
value, while reproducing the observed Higgs mass. According to the size and sign of the 
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one-loop contribution in (fT^ . the quartic coupling A must be smaller (larger) than the 
tree-level value for a positive (negative) loop contribution. 

The triple Higgs coupling is the third derivative of flT^ . which can be simplihed as 


^hhh — SAi; -|- 


V 




aiV 


m, 


2 

mf 


+ 31og^ 


327r2 ^ 

i 

Using flT5]) and flTT)) . the triple Higgs coupling flT5]) can be simplihed as 


(15) 


, _ 

'^hhh — - 

V 


1 + 


V 


96n‘^m\ 


E 


mj 


(16) 


The size of the radiative effects can be parameterized by the following dimensionless 
quantities: 


= 


fj? — Xv‘^ 


^ n m? — 2Xv‘^ 
= ^^—2 -: 


^Xhhh — 


Xhhh QXv 
Xh 


(17) 


/i rnii ^hhh 

These measure the relative strength of the radiative contributions to the Higgs bare mass- 
squared, the physical mass-squared jj?, and the triple Higgs coupling, respectively. Using 
the previously mentioned benchmark points, in Fig. |2]we plot the triple Higgs coupling 
versus the mass-squared parameter jj?, for both triplet and quintuplet models. We also 
show the relative strength of the radiative contributions to the Higgs mass, triple Higgs 
coupling and the parameter /i^, as dehned in Eq. flTTI) . 

One notices from Fig. [2]- Top that the mass-squared parameter /i^ can be large, even 
up to 100 (500) times the Higgs VEV-squared v'^ for the triplet (quintuplet) models. The 
larger values are required in order to balance the radiatively induced mass term in the 
Lagrangian, i.e., the second term on the left-hand side of firoj) . The radiative corrections 
can also be negative, depending on the value of the Higgs quartic coupling; i.e., for 
A > 0.08. Due to the fact that the extra helds in the quintuplet model are much heavier 
than those of the triplet model, their radiative contributions are larger and therefore the 
/i^ parameter values are larger, as it is evident from the hgures. From Fig. |2]-Bottom, 
one notices that the relative radiative-contributions to the Higgs mass and triple coupling 
are proportional, i.e. when the Higgs mass is completely generated radiatively, the triple 
Higgs coupling is also dominated by radiative effects. We also observe that for most of the 
benchmark points, in both the triplet and quintuplet models, the mass-squared parameter 
/i^ is fully radiative, as shown in the palette. 

The relevant quantity for collider phenomenology is the relative enhancement in the 
triple Higgs coupling, with respect to the SM value, which is dehned as 

Xhhh 


A = 


_ \SM 
'^hhh 
\SM 
'^hhh 


(18) 











Triplet Model 


Quintuplet Model 


0.16 




Figure 2: Top: the triple Higgs coupling versus the mass parameter parameter p?, both in 
units of the Higgs VEV, for the triplet (left) and guintuplet (right) models. The palette 
gives the Higgs guartic coupling and the red line shows the SM triple Higgs coupling 
value. Bottom: the relative radiative contribution to the triple Higgs coupling versus the 
relative radiative contribution to the Higgs mass. The palette gives the relative radiative 
contribution to the Higgs mass parameter . 


According to Eq. ffTBD . the relative enhancement of the triple Higgs coupling is given by 
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i^SM m; 


riia'r 

* I 

T" 


967r^r 




all 


riia'r 

-2” 


(19) 


In Fig. m we show the relative enhancement of the triple Higgs coupling, Eq. ([T9]), for the 
benchmark sets used previously. The hgure shows that the relative enhancement of the 
triple Higgs coupling, with respect to the SM, are larger for large values of the quartic 
coupling A, and smaller for small values of the charged scalar mass. Also, one notices that 
the relative enhancement in the triple Higgs coupling, with respect to the SM value, is 
always positive, contrary to other models 113120], and furthermore it can exceed 35% for 


9 














Triplet Model 
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Figure 3: The relative enhancement in the triple Higgs coupling with respect to the SM, 
A, versus the Higgs guartic coupling. The palette shows the mass of the charged scalar 
in GeV. The black point at (X = X^^,0) refers to the SM. 

both the triplet and quintuplet models. 

4 Electroweak Phase Transition 

The SM cannot successfully explain baryogenesis 121 ] for two reasons: (1) the CP violating 
source in the CKM matrix is too small and (2) the electroweak phase transition (EWPT) 
is not strongly hrst order. The latter is required to suppress the B+L violating processes 
in the broken phase, inside the bubble, when its wall is expanding during the transition. 
In the SM, the criterion for a strongly first-order EWPT 1221 . 

v{r)/T, > 1, (20) 

is not fulhlled since the ratio is given by Vc/T^ ~ A, which would require a Higgs mass 
below 42 GeV |23]. Here T^. is the critical temperature at which the effective potential 
exhibits two degenerate minima, one at zero and the other at v{Tc). Both Tc and v{Tc) 
are determined using the full effective potential at finite temperature, which is given by 


K/G + V.in,(.h,Ty, (21) 

i 

roo 

JB,F{a)= / log(l =F exp(-\/x2 -h a)). (22) 

Jo 

In the above, we include an important leading term from the higher-order loop corrections, 
which can play an important role during the EWPT dynamics, namely the so-called daisy 
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contributions [25] 




(23) 


Vring{h,T) = 


I 


The summation is over scalar and longitudinal gauge degrees of freedom, with fnf{h, T) = 
mj{h) + n(T) their thermal masses, and n(T) are the thermal parts of the self energy 
(given in the appendix). For our analysis we include the daisy contributions by following 
an alternative approach to Eq. fl25]) . i.e. by replacing the held dependent masses of the 
scalar and longitudinal gauge helds by their thermal masses in the full effective 

potential fl2T|) . In order to account for all the (heavy and light) degrees of freedom, we 
evaluate the integrals fl22|) numerically. 

The strength of the EWPT can be improved when new bosonic degrees of freedom 
are present, as occurs in the present models. It is clear from 0141) that for large values 
for the couplings {Xsh, Xht, Xht for triplet model and Xsh, Xhti,2 for quintuplet model} 
and/or small mass-values for the extra (singlet and multiplet) scalars, the one-loop correc¬ 
tion to the Higgs mass can be signihcant, allowing the Higgs self-coupling to be smaller. 
Consequently one can fulhll the criterion fl2U]) without conhicting with recent Higgs mass 
measurements m 

Analyses of similar models iminiEe] has shown that extra scalars can help to generate 
a strongly hrst order EWPT by: (a) relaxing the Higgs self-coupling A to be as small as 
(P(10“^); and (b) enhancing the value of the effective potential at the wrong vacuum at 
the critical temperature, without suppressing the ratio v{Tc)/Tc, which relaxes the severe 
bound on the mass of the SM Higgs. 
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Figure 4: The ratio v{T^/Tc versus the critical temperature for the (left) triplet and (right) 
quintuplet models. The palette shows the Higgs quartic coupling X. 
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In Fig. m we plot the ratio v{Tc}/Tc verses the critical temperature T^, using the 
20,000 benchmark points for the triplet and quintuplet models. The hgure shows that the 
EWPT is strongly hrst-order for a majority of the benchmark sets, with the ratio v{Tc)/Tc 
predominantly taking values between 1.2 and 1.5 in both the triplet and quintuplet models. 
The transition temperature is a bit larger than the typical SM value ~ 100 GeV, and can 
be as large as 150 GeV and 140 GeV for the triplet and quintuplet models, respectively. 
One can read from the palettes in Fig. 0] that, for hxed Higgs quartic coupling values, the 
ratio v{Tc)/Tc is inversely proportional to the critical temperature. 

Inspecting Fig. 01 one is lead to the conclusion that the increased EWPT strength is 
not only a consequence of a small Higgs quartic coupling A, but can also be due to the 
transition dynamics; the existence of heavy scalars makes the Higgs VEV slowly decaying 
with respect to the temperature. Gonsequently the evolving (increasing or decreasing) 
effective potential at the wrong vacuum makes the transition occurring at the mentioned 
temperature values, therefore giving a large ratio v{Tc)/Tc. 

An interesting issue, discussed in the literature 123. is a possible correlation between 
the EWPT strength and the relative enhancement in the triple Higgs coupling A. In 
Ref. [in] it was shown that such a correlation is not clear for a model with extra charged 
and neutral scalars from two inert doublets. In Fig. 0] we plot the relative enhancement 
of the triple Higgs coupling. A, versus the ratio v{Tc)/Tc, for the 20,000 benchmarks for 
the triplet and quintuplet models. From Figure |5l it is not clear whether a correlation 
between the relative enhancement in the triple Higgs coupling and the EWPT strength 
exists. This issue deserves a detailed and model-independent investigation. 

5 The Higgs decay channels /i — )► 77 and h 7Z 

In July 2012 , the ATLAS [28] and GMS [29] collaborations announced the observation of a 
scalar particle with mass ~ 125 GeV, with roughly 5a conhdence level. Subsequently this 
value was updated to rrih = 125.09 ± 0.21 GeV d- An important question is whether or 
not this really is the SM Higgs or an alternative Higgs-like state with different properties. 
Indeed, a ht to the data, performed by both the ATLAS and GMS collaborations, seems 
to (almost) show agreement with the SM, with the reported values being 1.17 =F 0.27 [30] 
and 1.13 ± 0.24 [21], respectively. 

Dehning as the branching ratio of h ^ 77 scaled by the SM value, we hnd that 
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Figure 5: The relative enhancement in the triple Higgs coupling A versus the ratio v{Tc)/Tc 
for the benchmark points used previously. 


the present models give 

_ B{h^ 77) 
"77 ^ 77) 


E, (r.) 

2 Ar M + in) 


(24) 


where i stands for all charged scalars Xi, Ti = with mxi being the mass of 

the charged particle Xi running inside the loop, iVc = 3 is the color number, and Qt is 
the top quark electric charge in units of |e|. The parameters 'di are given for triplet and 
quintuplet members in Table [T] In the above, the loop amplitudes Af"' for spin 0, spin 
1/2 and spin 1 particles are given by [32] 


i^) = -^ f (a;)], 

^V 2 (^) = [x + ix-l)f (x)], 

A'P (x) = —x~‘^ [2x‘^ + 3a; + 3 (2a; — 1) / (a;)] , (25) 


fix) 



arcsin^ i^/x) 

IH 

[logjtdis «] 


X < I 

X > 1. 


(26) 


Another important Higgs decay channel that can be modified by the presence of extra 
charged scalars is h —)■ 'jZ. This channel is similarly parameterized as 


R.yz — 


Bih^-fZ) 

BSMih ^ 


1 + 




'ilX (tw, Cw) + 2 (1 - 8s2/3) (n, (t) 


(27) 
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^W^-W 

T+ 

‘^^HT1+^HT2 

2 
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Cyj 2 

j^++ 

4:XhtI + ^HT2 

4AffTl+AffT2 

2 

j^+++ 

9Xhti 

6-9*™ \ 
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Table 1: The parameters ■(9* and Xj, which are relevant for the Higgs decay channels 
h —)■ 77 and h ■jZ. 


where (i = and the functions are given by [32] 


y) = h (a;, y ), 

^112 v) = v) - ^2 {x, y ), 

y) = [(1 + 2a;) si/cl - (5 + 2x)\ h {x, |/) + 4 (3 - sl/cl) h (x, y ), 

h {X, y) - -2^ + 2 (.-yf + (.-yf ’ ^2 {X, y) - ^(cc-y) ’ 


9[x) = 


's/I—' 


2(x-y) {x-y) 

\/x~^ — 1 arcsin {^/x) 


X < 1 

X > 1. 


(28) 


(29) 


The parameters x* are shown in Table [T] for both the triplet and quintuplet models. 

The deviation of the channels h —?■ 77 , 7 Z from their SM values is sensitive to the 
mass of the scalars and the strength with which they couple to the Higgs doublet, i.e. on 
the parameters m|, ml^, \sh, ^ht and \ht for the triplet model, and on m|, ml, 
ml^, ml^_^_, XsH and Xhti,2 for the quintuplet model. Depending on the relative sign 
of the couplings to the Higgs doublet, the new contributions can strengthen or weaken 
the deviation of B{h 77 ) from its SM value. In Fig. El we present versus R-yz for 
the considered 20,000 sets of benchmark parameters for both the triplet and quintuplet 
models. 

We remark that some benchmarks in the triplet model are already excluded by the 
recent measurements of ATLAS mi and CMS [21], while more precise measurements 
are required to probe benchmarks for the quintuplet model. In contrast to other models 
with extra singlets [1] and doublets [19], the decay channel h —)■ 7 Z can be significantly 
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Figure 6: The modified Higgs decay rates B{h ^ 77) vs B{h —?■ ■yZ), scaled by their SM 
values, due to the extra charged scalars, for 20,000 randomly chosen sets of benchmark 
parameters for the triplet (red) and guintuplet (green) models. The intervals between the 
magenta (green) lines represent the ATLAS (CMS) recent measurements on the /i —)■ 77 
channel, and the blue point represents the SM. 

modified, with respect to the SM value, particularly in the quintuplet model. This can be 
understood from the large Xj coupling values for the scalar multiplet members in Table [H 

6 Conclusion 

Models of radiative neutrino mass with DM candidates can explain some of the short¬ 
comings of the SM while generating observable experimental signals. In this work, we 
performed a detailed study of the scalar-sector phenomenology for a pair of three-loop 
neutrino mass models with DM candidates. The models, referred to as the triplet [6] and 
quintuplet [7] models, generate neutrino mass via a diagram with the same topology as the 
KNT model. We investigated the effect of the extra scalars on the Higgs mass, the triple 
Higgs coupling, and the Higgs decay channels h —)■ 77,7Z. We also studied the strength 
of the electroweak phase transition. In both models, it was shown that the beyond-SM 
multiplets can modify the triple Higgs coupling and the Higgs decay channels away from 
their SM values. The electroweak phase transition was found to be strongly hrst-order in 
signihcant regions of parameter space. Measurements of the Higgs decay channels already 
exclude some regions of parameter space for the triplet model, and future improvements 
will further explore the parameter space for both models. 
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A Field Dependent Masses 


I charged scalar and SM field-dependent masses are given by: 


The _ 

2 

2 2 I \ 7,2 I TT 2 2 I o \ 1,2 I TT 2 Vt i 2 2 2 92 

m^ = -jj + Xh mf^ = -n + 3Ah -h H^, = —h , = — 

2 \ 

2 _ •S'l 1,2 I TT 2 _ 929i 7,2 2 _ 2 I '^SH ,2 , tt 

^ + ns, 


— + n^v, 

(30) 


with 

hljs = (l2A + + 'ig'^ + 3?/^ + 2Xsh) —, 715 = (4A5_h- + 4 A 5 + 3g'^) —, 

1 1 ov 

(31) 


rri _ rri _ 11^20^2 ttL _ ^„/2m2 ttT _ ttT _ q 
~ ‘-^Ws ~ g y ) ‘-^B ~ igi/ ~ ^^B ~ 


Here, we ignored the triplet and quintuplet contributions since they decouple from the 
thermal plasma due to their large masses, relative to the relevant typical temperature of 
0(100 GeV). 

The triplet members field dependant masses are given by: 

XhT - XhT ,2 , TT 2 2 , ‘^^HT - Xht ,2 , tt 

- - h + Hr, m_^_ = /if H-^- h -h H^, 

—^h‘^ + n^, Hr = (9g^ + 3g'^ + 2 Xht + + ‘2Xst) (32) 


2 2 I 

rriQ = fij,+ 


m 


'++ 


— /^T + 


' ' ■ ^ 2 ' ' i. 

and the quintuplet members field dependant masses are given by: 

2 2 , ^hti + ^ht 2,2 , tt 2 2 , 4A/Tri + 3A/Tr2 7 2 , tt 

Tfl_ - jjjrp -\ --- h -\- li(P, ?7 Zq - jlrp -\ --- h 

2 o 

dA^n + ^HT2 7 2 , TT 

- - - h + Ht, 


2 2 , 

Tnj_ - jjjrp -\- 

.2 


2AhT1 + AhT2,2 , tt 2 _ 2 , 

- - - h Ht, ^++ — /^T T 


= ^^T + ^ 


4 

—-|- Ht, Ht = {9g^ + 3g'^ + 2 Xst + 4Ai5Ti + 2Xht2) 
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